function u = func_ex2(n, m, i, j)
    global bounds
    
    x = 0 : bounds(1) / (n-1) : bounds(1);
    y = 0 : bounds(2) / (m-1) : bounds(2);
    
    u = y(j)*(x(i)*(2 - x(i)))^(1/2);
end



